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ABSTRACT 

We performed a new and accurate fit of light and radial velocity curves of the Large Magellanic 
Cloud (LMC) Cepheid -OGLE-LMC-CEP-0227- belonging to a detached double- lined eclipsing binary 
system. We computed several sets of nonlinear, convective models covering a broad range in stellar 
mass, effective temperature and in chemical composition. The comparison between theory and obser- 
vations indicates that current theoretical framework accounts for luminosity -V and I band- and radial 
velocity variations over the entire pulsation cycle. Predicted pulsation mass -M = 4.14 ± O.OBMq- 
and mean effective temperature -re=6100±50 K- do agree with observed estimates with an accu- 
racy better than \a. The same outcome applies, on average, to the luminosity amplitudes and to 
the mean radius. We find that the best fit solution requires a chemical composition that is more 
metal-poor than typical LMC Cepheids (Z=0.004 vs 0.008) and slightly helium enhanced (Y=0.27 
vs 0.25), but the sensitivity to He abundance is quite limited. Finally, the best fit model reddening 
-E{V — I) = 0.171 ± 0.015 mag- and the true distance modulus corrected for the barycenter of the 
LMC -/io,LMC=18.50±0.02±0.10 (syst) mag-, agree quite well with similar estimates in the recent 
literature. 

Subject headings: stars: variables: Cepheids — binaries: eclipsing — Magellanic Clouds — stars: 
distances 



1. INTRODUCTION 

Classical Cepheids are among the most used primary 
distance indicators, currently at the base of an accu- 
rate extragalacti c distance scale (iFreedman et al.l [20011 : 
iSaha et al.l 120011: IFreedman et al.ll20lH) . thanks to their 
characteristic Period-Luminosity (PL) relation. At the 
same time, these intermediate mass (typically from 3 
to 13 Mq) helium burning pulsating stars are used 
as tracers of relatively young {^ 10 -;- 100 Myr) stel- 
lar populations. The physical basis of the PL relation 
is well known. By combining the period-density and 
the Stephan-Boltzmann relations, we obtain a Period- 
Luminosity-Color-Mass (PLCM) relation. Stellar evo- 
lution predicts the occurrence of a Mass Luminosity 
ML) relation fo r intermediate-mass He burning stars 
Bono et "aLll2000l ) and the PLCM can be easily converted 
into a Period-Luminosity-Color (PLC) relation holding 
for each individual Cepheid. The projection of the PLC 
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onto the PL plane gives the PL relation. On this basis, 
the investigation of Cepheid pulsation properties is also 
crucial for providing independent constraints on the ML 
relation. 

In this context one of the most popular open prob- 
len i is the s o -calle d mas s discrepancy problem suggested 
by iChristvl (|T970l ) and iStobid (|l969l ). The quoted au- 
thors noted that the Cepheid masses evaluated from 
the comparison of theoretical isochrones with obser- 
vations were larger by almost a factor of two than 
the masses obtai ned from applicati on of a period-mass- 
radius relation (jFricke et al.l Il971| ) . Even if a reduc- 
tion of this discrepancy between evolutionary and pul - 
sational masses was obtained by iMoskalik et al.l ()1992f ) 
on the basis of updated radiat ive op acities, namely 
by the OPAL (Iglcsias & Rogcri |199i| and the Opac- 
itv Project (fSea ton et a l..il994) g roups, several authors 
(iCaputo et air2 005: Na tale et al. 2008: B ono et al.l[200l 
'Kelle r fc Wo'odri2006t lEvans et al.l i2005tlBrocato et al.l 
2003) noted tha the effect remained of the order of 10%- 
20% for both Galactic and Magellanic Cepheids. 

A crucial input to this scenario was provided by 
the recent detection of two Cepheids in well de- 
t ached, double-lined eclipsing b inary systems in the LMC 
(Piet rzvhskieFaLl l2010l l20ll OGLE-LMC-CEP-0227, 
OGLE-LMC-CEP-1812). The geometry and the preci- 
sion of both photometric and spectroscopic data allowed 
the authors to measure the masses of these Cepheids with 
the unprecedented precision of 1%, thus fixing a corner- 
stone i n the solution of the ma ss discrepancy problem. 
Indeed, iCassisi fc Salarij ()2011[ ) were able to find evolu- 
tionary masses consistent with the dynamical estimate 
for OGLE-LMC-CEP-0227, on the basis of stellar evolu- 
tion models that include a moderate amo unt of convec- 
tive core overshooting. On the other hand, iNeilson et al.l 
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(|201lD found that both moderate convective core over- 
shooting and pulsation-driven mass foss are required 
to achieve a solution of the m ass discrepancy problem 
an dPrada Moroni et al.l ()2012l ) adopt a new Bayesian ap- 
proach to constrain the intrinsic stellar parameters of 
OGLE-LMC-CEP-0227, in particular its mass and age, 
by varying the efhciency of mass loss and core overshoot- 
ing. 

In this investigation we constrain the pulsation proper- 
ties of OGLE-LMC-CEP-0227, by modelling its light and 
radial velocity curves with nonlinear convective models 
specifically computed for this purpose. In §2 we present 
the observational and theoretical framework, while in §3 
we discuss the procedure adopted to fit light and velocity 
curves. In the same section, we also investigate the im- 
pact of metal and helium content on the best fit models. 
Conclusions and final remarks close the paper. 

2. OBSERVATIONAL AND THEORETICAL FRAMEWORK 

The observables we plan to use in the comparison be- 
tween theory and observations are the V and I-band light 
curves (with the amplitudes Ay, Aj) and the radial ve- 
locity curve. Note that we define the light and the ra- 
dial velocity curves observables even if they have been 
estimated by deblending the cumulative light and radial 
velocity curves of the binary system. To further con- 
strain the projection factor, i.e. the parameter adopted 
to transform radial velocity into pulsation velocity, the 
radial velocity amplitude in our approach was assumed as 
a free parameter. The above observables together with 
the pulsation periods are available for a large number 
of classical Cepheids, in particular for those observed 
with one of the different flavors of the Baade-Wesselink 
method. 

During the last few years it has been found that 
nonlinear, convective Cepheid models can play a cru- 
cial role ill fitting the quoted observab les (JBono et"al] 
[19991 [20001 [200a Keller fc WoodI [20061) . These mod- 
els, once fixed the ML relation, provide very accurate 
fits of the above observables. On the basis of this 
theoretical framework and by assuming plausible val- 
ues for the chemical composition (metal -Z- and he- 
lium -Y- abundances) and for the effective tempera- 
ture very accurate fits have been provided not only 
for Galactic and Magellanic Cloud (MC) Cepheids (see 
e.g. 'Wood et al."1997; 'Bono et al."2002"; Kell er fc Wo"odl 
[2006; Natalc et al. 2008; Marconi ct al. 2013,), but also 






for other groups of radial variables (iMarconil |2009| ) . A 
further advantage of this nonlinear pulsation approach is 
that it provides independent estimates of the individual 
distance modulus and of the reddening (Bono ct al."2002'; 
Marconi & Clcmentini 2005; Marconi & Dcgl'Innoccnti 
20071: iMarconi et al.ll201,'^ . 



Moreover, intrinsic parameters of radial variables 
based on pulsation observables can be compared with 
those based on evolutionary observables. However, this 
approach is partially hampered by the fact that pulsa- 
tion models are envelope models, and therefore, they re- 
quire the mass-luminosity relation predicted by evolu- 
tionary models. The binary LMC Cepheid OGLE-LMC- 
CEP-0227 provides, for the first time, the opportunity 
to break this degeneracy. The unprecedented accuracy 
of the dynamical mass allows us to use a restricted range 
in mass values. Moreover, we also have for the same 



object a precise estimate of the mean effective tempera- 
ture, this means that we can change the mean luminos- 
ity until predicted and observed period agree with each 
other. We still lack detailed information concerning the 
chemical composition, but for LMC Cepheids spectro- 
scopic investigations are already available in the litera- 
ture (Romaniello et al. 2008). The observed parameters 
of OGLE-LMC-CEP-0227 are hsted in the first row of 
Table 1 0. 

3. FIT OF LIGHT AND VELOCITY CURVES 

To perform an accurate fit of the pulsation observ- 
ables -V and I-band light and radial velocity curves, we 
constructed a large set of pulsation models by adopt- 
ing the typical chemical composition of L MC (Z=0.008 , 
Y=0.25) young stella r population (Luck et al.l 119981 : 
[Romaniello et al.ll2008l: iMucciarelh et al.ll201lD, th e dy- 
namical mass measured bv iPietrzvhski et al.l ( 2010[) and 
a broad range of effective temperatures with a step of 50 
K. The luminosity of individual models for each value of 
the effective temperature was changed in order to match 
the pulsation perioc0 (see Table 1). Fig. 1 shows the 
comparison between a subset of the computed models 
and the observations for OGLE-LMC-CEP-0227. Data 
plotted in this figure indicate that theory agrees quite 
well with observations. 

However, to constrain on a quantitative basis the qual- 
ity of the fit between the diiferent sets of models and 
observations, we adopted a x^ analysis. We phased the 
light and pulsational velocity curves of models in order to 
locate the maximum of the V-band light curve at phase 
zero. Moreover, the best fit between theory and obser- 
vations was computed by minimizing the following x^ 
functions: 
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From top to bottom they refer to light and radial velocity 
curves. Note that the index i runs over the two photometric 
bands and j over the number of measurements, Npoints= 68, 
420, 135 for V, I and radial velocity, respectively. The pre- 
dicted light and pulsational velocity curves are indicated with 
Mmod and Vmod and they have been evaluated at the same 
phase -4>j- of the f^ measurement -mj, Vj"'''- by fitting with 
a spline the predicted curves. The magnitude x^ function is 
minimized with respect to shifts in 5M and to the possible 
shifts in phase -5<j>- (typically ranging from 0.01 to 0.02 for 
the models listed in Table 1), while the velocity x^ function 
is minimized with respect to the 50 (the same as for the pho- 
tometry at the 0.01 level), the barycentric velocity -7- and 
the projection factor p, i.e. the parameter adopted to convert 
the observed radial velocity into pulsational velocity. The 

* The mean radius listed in column 11 is slightly different than 
the value given in Pietrzynski et al. (2010). The new estimate plus 
the projection factor -p- are based on the results by Pilecki et al. 
(2013, in preparation). 

" The exact match between observed and predicted period was 
fixed by using linear nonadiabatic pulsation models that provide 
the static envelope st ructure to the nonlinear hydrodynamical mod- 
els IjBono et al.ll999l ) . The period predicted by nonlinear models at 
limit cycle stability can differ from the linear one by a few percent. 
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Fig. 1. — From top to bottom V- and I-band light curves and radial velocity curve. Open circles mark observations, while red lines display 
nonlinear pulsation models constructed at fixed stellar mass (M=4.14M0) and chemical composition (Z=0.008, Y=0.25). From left to 
right the different panels show pulsation models constructed assuming different values of the mean effective temperature (labeled values). 
The luminosity of individual models was changed until predicted and observed period agreed with each other. The blue open circles mark 
the best fit solution, according to the total rms, for this set of models. The apparent magnitudes of the light curves were transformed into 
absolute magnitude using the true distance modulus and the reddening listed in Table 1. 
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phase difference applied to light and radiative velocity curves 
are quite similar, and indeed we found that both the observed 
and the predicted phase lag, defined by using Fourier phases, 
between V-band light and radial velocity curves is -0.2. It 
is worth noting that the shifts in SM provide the Cepheid 
true distance modulus and the reddening, while 7 gives the 
bary centric velocity of the system. 

We neglected the models for which the x analysis provided 
a shift in radial velocity -7- which are not consistent, within 
the errors, with the observed mean barycentric velocity (see 
Table 1). The rms for the different observables, given by the 
X^ analysis, have been normalized to their pulsation ampli- 
tudes and summed in quadrature to obtain a total rms. 

The best fit solution, for each set of models, satisfy the 
above requirements concerning the shifts and attains the 
smallest total rms. The best fit model is plotted with open 
blue circles in Fig. 1 and the intrinsic parameters are listed 
in Tab[T](see discussion below). 

By applying the x^ analysis to the models computed at 
fixed chemical composition (Z=0.008, Y=0.25) and stellar 
mass {M — 4:.14Mq), we find that the effective temperature 
of the best fit model -Te = 6050±50 K (see Fig. 1)- agrees 
quite well with the empirical estimate (Te = 5900±250 K). 
The same outcome applies to the other observables, since the 
agreement between theory and observations is better than la 
for the mean radius and the V-band amplitude, and of the 
order of 2cr for the I-band. Note that the uncertainties on the 
intrinsic parameters of the best fit model were estimated ac- 
cording to the step in mass value and in effective temperature 
of the different sets of models. 

To further improve the accuracy of the intrinsic parame- 
ters, we computed a new set of models at fixed chemical com- 
position and effective temperature (T^ = 6050 K). The new 
models were constructed assuming a step in mass of ~0.Q6Mq 
and once again the luminosity was changed to fit the observed 
period. The x^ analysis indicates that the best fit is given by 
a model with a stellar mass of 4.08Mq and a luminosity level 
logL/L0=3.13 dex. Data plotted in Fig. 2, where a subset 
of the computed models is shown, display that theory and 
observation agree quite well with each other. 

3.1. Dependence on chemical composition 

In order to constrain the dependence of the best fit so- 
lution on the adopted chemical composition, we computed 
additional pulsation models by varying either the metallicity 
or the helium content. In particular, we computed new sets 
of models at fixed helium content -Y=0.25- and three differ- 
ent metal abundances: Z=0.006, Z=0.004 and Z=0.002. For 
each set of models, we performed the same fit to the observed 
curves both at fixed mass and at fixed effective temperature, 
following the same approach we adopted for the models at 
canonical chemical composition (Z=0.008). The best fit solu- 
tions for Z ranging from 0.004 to 0.008 are plotted in Fig. 3 
and the corresponding stellar parameters for all the metallic- 
ities are listed in Table [l] 

The total rms of the best fit models for the four dif- 
ferent metallicities, indicate that the global best fit is 
given by the model with Z=0.004, Y=0.25, M=4.26±0.06 
Mq, log L/L© =3. 15±0. 02 and an effective temperature of 
6050±50 K, again in ig ood agreement with the results by 
IPietrzvnski et all (|2O10i ) (see also Table 1)13. Current analy- 
sis further confirms the sensitivity of pulsation models to the 
metallicity, and in turn the opportunity to constrain this pa- 

^^ The match of the model at Z=0.002 with the data predicts 
a b arycentric velocity -7— at odds with the mean value given by 
IPietrzvnski et al. (201(|) and for this reason it was not taken into 
account (see §3). However, it was kept in the analysis only for the 
sake of completeness. 



rameter with the comparison between predicted and observed 
light and radial velocity variations (Natalc ct al. 2008). 

To further constrain the impact that chemical composition 
has on the comparison between predicted and empirical ob- 
servables, we also computed new sets of pulsation models at 
fixed metallicity -Z=0.004- and three enhanced helium con- 
tents, namely Y=0.26, 0.27, 0.28. The dependence of the 
best fit solution on the helium content is shown in Fig. 4. 
Data plotted in this figure and predicted observables listed 
in Table [1] indicate that the best fit solutions are those with 
He enhanced (Y=0.27, Y=0.28) composition. However, the 
rms values clearly show that pulsation observables are less 
sensitive to changes in He than in metal content. On the ba- 
sis of this evidence and the fact that the set of models with 
the highest He content (Y=0.28) give bar ycentric velocity - 
7- sli ghtly larger than the value given by IPietrzvnski et aP 
l|2010D . we assume that our global best fit model is the one 
with Z=0.004 and Y=0.27. 

4. CONCLUSIONS AND FINAL REMARKS 

We computed new sets of nonlinear convective pulsation 
models in order to match the observed V- and I-band light 
curve, the radial velocity curve and the radius curve of the 
LMC Cepheid OGLE-LMC-CEP-0227. This variable be- 
longs to a detached, double-lined eclipsing binary system, 
and for the first time its dynam ical mass was es t imate d 
with an accuracy better than 1% IPietrzviiski et ahl ()2010l ). 
By adopting the canonical chemical composition for LMC 
Cepheids (Z=0.008, Y=0.25), the best fit model gives a 
value of the Cepheid mass -M = 4.08 ± 0.06Mq- and of 
the effective temperature -Te=6050±50 K- that agree quite 
well with the observed values. This outcome also applies 
to the mean radius (the difference is smaller than la) and 
to the pulsation amplitudes (the difference is smaller than 
2a). The true distance modulus corresponding to this so- 
lution is ^o=18.44±0.02 where the small uncertainty is the 
only contribution of the model fitting technique. By adding 
a conservative systematic error of 0.1 mag to account for 
several residual uncertainties introduced by the adopted tur- 
bulent convective model, the static model atmosphere, the 
input physics, together with observational errors, we found 
that current distance is in remarkable agreement with sev- 
eral independent estimates of the LMC distance in the recent 
literature (j Bono ot al. 2002; Marconi & Clcmentini 2003; 
Keller fc Wood 2 006: Testa ot al. 2007; Mohn aro et al. 20l3 : 
Marconi et al.l 12013: Piotrzviiski ct al. 2013). To provide a 



firm estimate of the LMC distance wc ev aluated the geomtri- 
cal co rrection by adopting the model by Ivan der Marel et al.l 
l|2002f) and we found that it is 0.062 mag. Therefore, cur- 
rent distance modulus once corrected for the barycenter of 
the LMC becomes /io,z,Afc=18.50±0.02 mag, (d=50.1±0.5 
kpc) that is in excellent agr eement with the very precise 
distance recently provided bv IPietrzviiski et al.l l|2013D using 
eight long-period, late-type double eclipsing binary systems 
(18.49±0.05 mag). The same agreement applies to the dis- 
tance modulus provided by iPrada Moroni et al. (2012) using 
evolutionary models (18.53±0.02 mag). 

However, we still lack an accurate measurement of the iron 
abundance of this Cepheid. Therefore, we also investigated 
possible variations of the chemical compositions. We found 
that the best fit model providing pulsation observables that 
agree quite well with observations and an even smaller rms 
is the model with Z=0.004, Y=0.27, M = 4.14 ± O.O6M0, 
log L/L0= 3.16±0.02 dex, Te^OlOOiSO K and i? = 34.3 ± 
0.2i?Q. This model provides a distance modulus corrected for 
the barycenter of the LMC of /io,lmc = 18.50±0.02 mag. 

Current pulsation mass estimate agrees quite well with the 
dynamical mass measurements. A similar agreement has also 
been found concerning the evolutionary masses suggesting 
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TABLE 1 

The first line give the parameters derived from the dynamical analysis BY IPiETRZYNSKI ET AL.I POIOI ') AND THOSE OF THE 

OBSERVED PHOTOMETRIC AND RADIAL VELOCITY CURVES, WHILE THE BEST FIT MODEL PARAMETERS ARE LISTED IN THE FOLLOWING SEVEN 

LINES OF THE TABLE: MeTALLICITY (COLUMN 1), HELIUM CONTENT (COLUMN 2), MASS IN SOLAR, UNITS (COLUMN 3), LUMINOSITY' (COLUMN 

4), Effective Temperature (column 5), Period (derived from the linear and the nonlinear analysis) in days (column 6), 
Amplitudes of the photometric curves expressed in magnitudes (columns 7, 8), V-I Reddening value in magnitudes (column 

9), DEREDDENED DISTANCE MODULUS IN MAGNITUDES (COLUMN 10), THE MEAN RADIUS (COLUMN 11), PROJECTION FACTOR (COLUMN 12), 

Barycentric Velocity in km/s (column 13) and Total RMS (column 14). 
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(mag) 


R 
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7 
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4.14 
±0.05 




5900 

±250 


3.793 


1.06 
±0.02 


0.63 
± 0.01 






33.7 
±1.5 


1.16 
±0.14 


256.70 
±0.10 




0.008 


0.25 


4.14 
± 0.06 


3.13 
±0.02 


6050 
±50 


3.793-3.805 


1.11 
±0.06 


0.69 
±0.04 


0.154 
±0.014 


18.45 
±0.03 


33.9 
±0.3 


1.08 
±0.05 


256.80 
±0.10 


0.077 


0.008 


0.25 


4.08 
±0.06 


3.13 
±0.02 


6050 
±50 


3.793-3.803 


1.12 
±0.09 


0.69 
±0.05 


0.154 
±0.012 


18.44 
±0.02 


33.7 
±0.2 


1.08 
±0.03 


256.78 
±0.10 


0.077 


0.006 


0.25 


4.14 
±0.06 


3.14 
±0.02 


6050 
±50 


3.793-3.818 


1.11 
±0.06 


0.69 
±0.03 


0.150 
±0.040 


18.47 
±0.02 


34.2 
±0.2 


1.15 
±0.03 


256.75 
±0.10 


0.074 


0.004 


0.25 


4.26 
±0.06 


3.15 
±0.02 


6050 
±50 


3.793-3.812 


1.11 
±0.05 


0.70 
±0.03 


0.150 
±0.020 


18.49 
±0.02 


34.5 
±0.2 


1.27 
±0.04 


256.79 
±0.10 


0.070 


0.002 


0.25 


4.02 
±0.06 


3.13 

±0.02 


6000 
±50 


3.793-3.812 


1.06 
±0.10 


0.68 
±0.05 


0.139 
±0.014 


18.46 
±0.03 


34.3 
±0.4 


1.39 
±0.01 


256.81 
±0.10 


0.078 


0.004 


0.26 


4.14 
±0.06 


3.14 
±0.02 


6050 
±50 


3.793-3.816 


1.10 
±0.06 


0.69 
±0.03 


0.156 
±0.013 


18.46 
±0.01 


34.4 
±0.2 


1.27 
±0.07 


256.78 
±0.10 


0.069 


0.004 


0.27 


4.14 
±0.06 


3.16 
±0.02 


6100 

±50 


3.793-3.805 


1.10 
±0.02 


0.68 
±0.01 


0.171 
±0.015 


18.44 
±0.02 


34.3 
±0.2 


1.20 
±0.08 


256.77 
±0.10 


0.068 


0.004 


0.28 


4.14 
±0.06 


3.15 
±0.02 


6100 
±50 


3.793-3.795 


1.08 
±0.02 


0.67 
±0.01 


0.172 
±0.015 


18.43 
±0.02 


34.2 
±0.2 


1.20 
±0.08 


256.81 
±0.10 


0.067 



either an increase in tiie efficiency of con vective core over- 
shooting during liydrogcn burning phases l|Cassisi fc Salari"i 
[2011; Prada Moroni ct al. 2012) and a pulsation driven mass 
loss inside the Cepheid instability strip (Neilson et al. 2011). 
Current pulsation and evolutionary results are hampered by 
the fact that they are only based on a single object and are 
affected by degeneracy among the different input parameters. 
The use of independent solid observables (period derivative, 
CNO abundances) will provide firm constraints on the actual 
evolutionary and pulsation status of the system (Neilson et 
al. 2012; Matthews et al. 2012). 

The p-factor given by our best fit solution -1.20±0.08- 
is consistent with the result obtained for the prototype 5 
Ce phei by[M crand et a.1.1 (|2005D and with the value suggested 



epr 

by JGroenewegenl (|2013l ) by adopting Galactic Cepheids with 
individual distances based on trigonometric parallaxes. The 
current estimate is also in very good agreement with the new 
binary best fit solution provided by Pietrzynski et al. (2013) 
and confirms the value predicted by Nardctto ct al. ( 2004) on 
the basis of Cepheid hydrodynamical models. However, the p- 
factor we found is smaller than the values recently determine d 
for short-period Cepheids in the LMC bv lStorm et all (|2011l '). 
The p-factor and its possibl e dependence on the pulsation pe- 
riod is still lively debat ed IjGieren et al.ll2005l : iGroeneweeenl 
120131 : iNgeow et alll2012l ) and we need to apply the same anal- 
ysis to several Cepheids before we can rich firm conclusions. 

The reddening estimate is also in very good agreement 
with recent evaluations based on reddening maps (Schlegel 
et al. 1998) and on evolutionary models (Prada Moroni et 
al. 2012). The estimated luminosity level is in remarkable 
agreement not on l y with the estimate of the binary modelling 
IPietrzviiski et akl (|2010D , but also with evolutionary prescrip- 
tions based on models accounting for mild convective core 
overshooting during central hydrogen burning and canonical 
mass loss rates (see figure 3 in Prada Moroni et al. 2012). 
This evidence indicates a very promising agreement between 
pulsation and evolutionary observables. Finally, it is worth 
mentioning that the nonlinear pulsation approach can provide 
firm constraints not only on the pulsation observables -at a 
few percent level- but also on the metal abundance. However, 



their sensitivity to the helium content is quite limited. 



We acknowledge financial support from PRIN INAF 2011 
(P.I. M. Marconi) and from PRIN MIUR 2011 (P.I. F. Mat- 
teucci). G.B. thanks ESO for support as science visitor. WG 
and GP thank for support from the BASAL Centre de As- 
trofisica y Tecnologias Afines (CATA) PFB-06/2007. Support 
from Foundation for Polish Science (program TEAM) and the 
Polish National Science Centre (program MAESTRO) is also 
acknowledged. RS is supported by Polish NCN (grant DEC- 
2012/05/B/ST9/03932). We dedicate this manuscript to the 
memory of R. Buchler. His suggestions and ideas in the field 
of stellar pulsation will be greatly missed. 



Marconi, et al. 






e -3 



-2.5 



I I I I I I I I 1 I M I I 1 I I 

4.02 M, 



(0 



B -3.5 



X 



> 




0.5 1 1.5 0.5 1 1.5 0.5 1 1.5 
Phase Phase Phase 



Fig. 2. — Same as Fig. 1, but the comparison between theory and observations was performed at fixed chemical composition (Z=0.008, 
Y=0.25) and effective temperature (Te=6050 K). From left to right the different panels display pulsation models constructed by assuming 
different values of the stellar mass (see labeled values). The blue circles show the best fit solution for this set of models. 
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Fig. 3. — Same as Fig. 1, but for pulsation models constructed by assuming the same helium content (Y=0.25) and different values of 
the metal abundance (see labeled values). The stellar mass and the effective temperature of the best fit solution, for each set of models, 
are also labeled. Blue circles display the best fit solution of the current set of models. 
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Fig. 4. — Same as Fig. 3, but for pulsation models constructed by assuming the same metal content (Z=0.004) and different helium 
abundances (see labeled values). The helium enhanced (Y>0.25) pulsation models cover a broad range in stellar mass and in effective 
temperature. The stellar mass and the effective temperature of the best fit solution, for each set of models, are also labeled. Blue circles 
display the global best fit solution. 
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